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Recent progress in cosmic ray physics covering the energy range from about 10^'' eV to 10^^ eV is reviewed. 
The most prominent features of the energy spectrum are the so called 'knee' at _E ~ 3 ■ 10^^ eV and the 'ankle' at 
few 10^** eV. Generally, the origin of the knee is understood as marking the limiting energy of galactic accelerators 
and/or the onset of increasing outflow of particles from the galaxy while the ankle is considered to mark the 
transition from galactic to extragalactic cosmic rays. Alternative theories do exist and shall be sketched. A key 
observable to answer the still open questions about the cosmic ray origin and to discriminate between various 
models is given by measuring the chemical composition or - more directly - by measuring energy spectra of 
individual cosmic ray mass groups. The status of present analyses is critically discussed and new experimental 
endeavors carried out in order to improve both the statistics and the quality of data particularly at energies above 
the knee will be summarized. 



1. INTRODUCTION 

The cosmic ray (CR) energy spectrum extends 
from a few hundreds MeV to above 10^" eV. 
Over this wide range of energies the intensity 
drops by more than 30 orders of magnitude. 
Despite the enormous dynamic range covered, 
the spectrum appears rather structureless and 
can be well approximated by broken power-laws 
dN/dE oc E"'. Up to energies of a few 10^'' eV 
the flux of particles is sufficiently high to en- 
able measurements of their elemental distribu- 
tions by high flying balloon- or satellite-borne ex- 
periments. Such studies have provided important 
information about the origin and transport prop- 
erties of CRs in the interstellar medium. Two 
prominent examples are ratios of secondary to 
primary elements, such as the B/C-ratio, which 
are used to extract the average amount of matter 
CR-particles have traversed from their sources to 
the solar system, and are relative abundances of 
radioactive isotopes, such as ^"Be to stable ^Be 
or ^^Al to stable ^^Al, which carry information 
about the average 'age' of CRs. With many new 
complex experiments taking data or starting up 
in the near future and with a possibly new gener- 
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ation of long flying balloons, this remains a vital 
field of research. 

Above a few times 10^''^ eV the flux drops to 
only one particle per square metre per year. This 
excludes any type of 'direct observation' even in 
the near future, at least if high statistics is re- 
quired. On the other hand, this energy is large 
enough so that secondary particle cascades pro- 
duced in the atmosphere penetrate with a foot- 
print large enough to be detected by an array 
of detectors on the ground. Such an extensive 
air shower (EAS) array typically has dimensions 
of a fraction of a square kilometre to more than 
1000 square kilometres and can be operated for 
many years to detect fluxes down to 1 particle per 
square kilometre per century or less. 

The most prominent features of the CR energy 
spectrum fall into the energy range covered by 
EAS experiments. The steepening of the slope 
from 7 9^ 2.7 to 7 ^ 3.1 at about 3 • lO^^ eV 
is known as the 'knee'. It was flrst deduced from 
observations of the shower size spectrum made by 
Kulikov and Khristianson et al. in 1956 [T] but 
it still remains unclear as to what is the cause of 
this spectral steepening and even as to what are 
the sources of the high energy CRs at all. At an 
energy above 10^* eV the spectrum flattens again 
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Figure 1. The all-particle CR energy spectrum weighted by showing the knee at 3 • 10^^ eV, a possible 
second knee at 10^^ eV, the ankle at about 3 • 10"'^^ eV, and the GZK-rcgion near 6 • 10^^ eV. References 
are given in the text. 



at what is called the 'ankle'. Because of the large 
size and/or magnetic field required to accelerate 
and confine charged particles above 10^* eV, the 
origin of CRs above the ankle is generally consid- 
ered to be of extragalactic (EG) nature. Finally, 
the question whether the spectrum extends be- 
yond the Greisen-Zatsepin-Kuzmin threshold of 
6 • 10^^ eV [2] is currently among the foremost 
questions in astro-particle physics as is reficctcd 
also by the number of presentations given at this 
conference. 

The main purpose of this paper is to review the 
experimental data in the energy range below the 
GZK-threshold, i.e. from about lO^^ to lO^^ eV. 
We shall discuss the energy spectrum, chemical 
composition, and anisotropics in their arrival di- 
rections and critically examine the astrophysical 
implications by taking into account the system- 
atical uncertainties of the data. 

2. THE KNEE REGION 

Mainly for reasons of the required power the 
dominant acceleration sites of CRs are generally 



believed to be shocks associated with supernova 
remnants (SNR). Naturally, this leads to a power 
law spectrum as is observed experimentally. De- 
tailed examination suggests that this process is 
limited to Eq/Z - 10^'^ eV [SU for standard 
galactic SNRs. This value can be extended up- 
ward with a number of mechanisms, for exam- 
ple by introducing higher magnetic fields, larger 
sources, quasi-perpendicular shocks, reaccelera- 
tion by multiple sources, etc. However, these 
assumptions and their effects are not free of de- 
bate and possibly, something more fundamental 
may be incorrect with the suggested supernova 
(SN) picture and its shock value Eq. In any case, 
if there is a typical maximum energy which de- 
pends linearly on Z for reasons of magnetic con- 
finement, then the spectrum of CR nuclei must 
become heavier with increasing energy as the hy- 
drogen cuts off first and then increasingly heavier 
nuclei reach their acceleration (or confinement) 
limits. 

A change in the CR propagation with decreas- 
ing galactic containment at higher energies has 
also been considered. This increasing leakage re- 
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suits in a steepening of the CR energy spectrum 
and again would lead to a similar scaling with the 
rigidity of particles, but would in addition predict 
anisotropics in the arrival directions of CRs with 
respect to the galactic plane. 

Besides such kind of 'conventional' source and 
propagation models (5|6j several other hypothe- 
ses have been discussed in the recent literature. 
These include the astrophysically motivated sin- 
gle source model of Erlykin and Wolfendale ^ 
trying to explain possible structures around the 
knee by a single recent and nearby SN, as well as 
several particle physics motivated scenarios try- 
ing to explain the knee due to different kinds 
of CR-interactions, e.g. by photodisintegration at 
the source '8' or by sudden changes in the char- 
acter of high-energy hadronic interactions during 
the development of EAS [9] . 

Recently, the 'Cannonball' model of CRs has 
been suggested as a radically different theory of 
CR origin [10]. It is inspired by mounting obser- 
vational evidence that, in addition to the ejection 
of a non-relativistic spherical shell, the explosion 
of core-collapse SNae results in the emission of 
highly relativistic bipolar jets of plasmoids of or- 
dinary matter, the 'Cannonballs' (CB). As the 
CB with a typical half of the Mercury mass propa- 
gates at relativistic speed through the interstellar 
medium, it encounters electrons, protons, and nu- 
clei kicking them up to higher energies elastically 
by magnetic deflection. These newly born CRs 
are then subject to propagation effects, similarly 
as in 'classical' theories. It is argued that this very 
simple concept explains all observed properties 
of non-solar CRs at all observed energies. There 
are two important differences to the conventional 
models: a) because of the specific kinematics of 
particle acceleration, the maximum energy of CRs 
(and thereby the knee positions) scale with the 
mass A of CRs rather than with their charge Z, 
b) since the CBs propagate rapidly from the inner 
SN and GRB realm of the Galaxy into its halo or 
beyond converting ISM particles to high energy 
CRs all along their trajectories, there is a much 
lower level of CR-anisotropy expected than in the 
traditional SN picture of CRs. 

Indeed, the low level of CR anisotropy even 
at energies above the knee is a long standing 



problem [TT|. Generally, the observed spectrum 
(piE) and the source spectrum Q{E) are consid- 
ered to be connected by a relation of the form 
(f){E) = Q{E) X TesciE). A simple power-law 
fit of the escape time to the available data gives 
Tcsc{E) oc E^^ with 6 « 0.6. Extrapolating Tcsc 
to 10^^ eV, for example, would lead to a value al- 
most as small as the light travel time across the 
galactic disk, implying a much larger anisotropy 
than is observed. 

From the discussion above it is obvious, that 
an answer to the question about the origin of the 
knee is of key importance to reveal information 
about the origin of galactic cosmic rays in gen- 
eral. Experimental access to such questions is 
provided by measurements of charged cosmic rays 
(the classical nucleonic component) and 7-rays 
by experiments above the atmosphere, and by 
the observation of air showers initiated by high- 
energy particles in the atmosphere. 

A wealth of information on potential cosmic- 
ray sources is provided by recent measurements 
of TeV 7-rays employing imaging atmospheric 
Cherenkov telescopes, most notably from the 
H.E.S.S. experiment. Their observation of the 
morphologies and energy spectra of the shell type 
SNRs RX J1713.723946 pTT3] and RX J0852.0- 
4622 14J are well in agreement with the idea of 
particle acceleration in the shock front. The spec- 
tra extend up to energies of 10 TeV and provide 
evidence for the existence of particles with en- 
ergies beyond 100 TeV at the shock front that 
emerged from the supernova explosions. How- 
ever, an unequivocally proof for acceleration of 
hadrons is still missing and questions arise also 
about the low number of established SNRs show- 
ing TeV 7-ray emission. For example, a recent 
Galactic plane survey of H.E.S.S. [15] reveals no 
SNRs brighter than these two in the region cov- 
ered. This apparent deficit of TeV-bright SNRs 
may pose some problems in explaining the high 
energy budget of galactic CRs. Remember that 
about 10 % of the mechanical energy released by 
the population of Galactic supernovae needs to 
be converted into CRs if all SNRs are sites of 
CR acceleration. Any reduction in the number of 
TeV-bright SNRs needs to be compensated for by 
a corresponding factor in the already large value 
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of the CR acceleration efficiency. 

To undoubtedly establish SNRs as the sites of 
CR acceleration and in order to constrain the con- 
ventional SN acceleration model from other pro- 
posed mechanisms, precise measurements of the 
primary CR energy spectrum and particularly of 
the mass composition as a function of energy are 
needed. Significant progress has been made here 
as well in recent years, but the situation is far 
from being clear. 

2.1. Comparison of direct and indirect 
measurements 

Cosmic ray measurements on balloons and 
spacecraft have an important advantage over 
ground-based air shower experiments: They de- 
tect the primary CR particles and measure its 
charge directly. This is because spacecraft experi- 
ments perform the measurement above the atmo- 
sphere and balloon-borne experiments typically 
perform their measurements with residual atmo- 
spheres of only ~ 5-10 g/cm^. This is a relatively 
small value compared to the typical hadronic in- 
teraction length of A/ ^ 90g/cm^ so that correc- 
tions for interactions above the instruments are 
of minor importance, at least for light particles, 
such as protons and He nuclei. This advantage is 
paid for at the expense of lacking statistics at high 
energies. For example, the largest of the current 
generation of balloon-borne detectors, TRACER 
[16j . reaches a sensitive volume of 2 x 2 x 1.2 m'^. 
It has been flown successfully for 14 days expo- 
sure from the Antarctic in 2003 and from Swe- 
den in summer 2006. The first 14 days flight 
time resulted in an exposure of ~ 75 m^ sr days 
and allowed to measure e.g. oxygen nuclei up to 
- 320 TeV and iron nuclei up to ~ 70 TeV. 

The largest exposure of all direct experiments 
has been reached by the Japanese American Co- 
operative Emulsion Experiment JACEE [TT and 
the RUssian-Nippon JOint Balloon collaboration 
RUN JOB [IH]. JACEE flew a series of thin 
('^ 8.5 radiation lengths) emulsion/X-ray film 
calorimeters on 15 flights during 1979-1994 and 
has reached an exposure of ~ 664 m^ hrs from 11 
analysed flights. Taking the zenith angle accep- 
tance out to ta,nd ^ 72—79° into account, this re- 
lates to approx. ^ 80 m^ sr days. RUNJOB flew 
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Figure 2. Proton (top), helium (centre), and iron 
(bottom) spectra from direct experiments com- 
pared to EAS data (based and updated from |T^). 
The single diagonal error bar in each panel in- 
dicates the eflFect of a ±15% uncertainty in the 
energy scale. 



roughly a similar set of X-ray films and emul- 
sion chambers on a series of 10 successful balloon 
fiights during 1995-1999 with a total exposure of 
575 m^ hrs. Both experiments were able to recon- 
struct proton spectra up to almost 1 PeV. 

Figure[2] shows a collection of the proton, he- 
lium, and iron spectra obtained by various ballon- 
and satellite-borne experiments compared to data 
from ground based experiments. Obviously, data 
from direct experiments are sparse above 100 TeV 
and uncertainties become very large with increas- 
ing energies, particularly for primaries heavier 
than protons. Reasonably good agreement be- 
tween RUNJOB and JACEE is observed in case 
of the proton spectrum, but the He-flux mea- 
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sured by RUNJOB is about a factor of two lower 
compared to other experiments. Comparing the 
slopes of the p and He spectra yields power law 
indices of ~ 2.7-2.8 for both elements in the en- 
ergy range 10-500 TeV/nucleon. The iron spec- 
trum appears somewhat flatter, 7fc — 2.6, par- 
ticularly when taking into account the extrapo- 
lation to the EAS data. Such a dependence of 
7 could be explained by charge dependent effects 
in the acceleration or propagation process. For 
example, non- linear models of Fermi acceleration 
in supernovae remnants predict a more efficient 
acceleration for elements with a large AjZ ratio. 
However, uncertainties may still be too large to 
allow for definite conclusions about differences in 
the acceleration and propagation mechanisms of 
different primaries. For illustration, the effect of 
an assumed (and possibly underestimated) uncer- 
tainty of ±15% in the energy scale is shown by 
the single error bar in each of the panels. 

It is remarkable to see direct measurements and 
EAS data starting to overlap each other. Clearly, 
EAS data below about 10^^ eV are dominated 
by systematic uncertainties while direct measure- 
ments suffer from statistical ones. With these 
caveats kept in mind, the agreement is very good. 
The EAS data of KASCADE [20], also shown in 
figure [21 have been reconstructed based on two 
different hadronic interaction models employed in 
the EAS simulations. Except perhaps for iron, 
the uncertainties caused by the interaction model 
are of similar size or even smaller than systematic 
uncertainties between experiments like JACEE 
and RUNJOB. Also shown in figure [2| are proton 
and helium spectra derived from emulsion cham- 
bers and burst detectors operated within the Ti- 
bet II air-shower array [21]. The results are in 
rough agreement with the KASCADE data. For 
reasons of clarity, only the results based on simu- 
lations with the CORSIKA [22] / QGSJet-model 
[53] are included for the Tibet data. Those ob- 
tained based on Sibyll [23] are similar within their 
error bars. There are some important peculiari- 
ties of the Tibet AS7 analysis to be pointed out 
here. The data are compared to EAS simula- 
tions assuming in one case a heavy dominated 
(HD) and in another case a proton dominated 
(PD) composition. In the HD-model a rigidity 



dependent knee E\ — Z ^ 1.5 • lO^'' eV is assumed 
and in the PD-model all mass components are as- 
sumed to break off at i?k = 1-5 • 10^^ eV. These 
assumptions are surprising, since no experiment 
ever has observed at break in the spectrum at 
such low energies. Furthermore, the experimen- 
tal data of Ref. [21] start only at energies above 

~ 4 • lO^^'eV, i.e. significantly above the as- 
sumed knee position. Moreover, because of in- 
sufficient separation power between proton and 
helium primaries, the authors have deduced the 
proton spectrum first by using a neural network 
algorithm. Next, the proton + helium spectrum 
has been reconstructed and, finally, the helium 
spectrum has been obtained by subtracting the 
number of proton events obtained in the first task 
from the proton + helium dataset obtained in the 
second task. Clearly, there are huge correlated 
errors to be expected in the helium spectrum de- 
duced that way. Also, it is not clear how the 
results depend on the ad-hoc assumptions made 
about the knee position. Because of the steeply 
reconstructed proton and helium spectrum, the 
authors then conclude, that the main component 
responsible for the change of the power index of 
the all-particle spectrum around 3-10"'^^ eV is com- 
posed of heavy primaries. However, there is no 
proof to this statement as the experiment is al- 
most blind to heavy particles (detection efficiency 
of iron « 4 %) . 

To conclude this topic, despite some contro- 
versy a reasonably good agreement between di- 
rect and EAS experiments has been achieved in 
recent years. At present, EAS experiments at 
their threshold energies are limited purely by sys- 
tematic uncertainties, while direct measurements 
suffer mostly from lacking statistics but also from 
systematic uncertainties in determining the ab- 
solute energy scale. There is some hope that 
new EAS experiments located at very high alti- 
tude will be able to push the measurements down 
to lower energies and at the same time also re- 
duce their systematic uncertainties. Direct ex- 
periments, on the other hand, may be able to in- 
crease their exposure at high energies. However, 
given the very steeply falling spectrum, it appears 
unlikely that balloon experiments will be able to 
extend the range of measurements beyond 1 PeV 
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any time in the near future. Thus, the chance of 
detecting the knee with direct measurements of 
protons to iron on baUoons is not hkely to occur 
without significant increases in the payload and 
flight duration capabihties of high altitude bal- 
loons. Even with 50 times the present JACEE 
p-He exposure one would still be unable to make 
definitive measurements about a break in the en- 
ergy spectrum beyond 200 - 300 TeV [25]. 

2.2. Air shower data at the knee 

As can be seen from figure [U a wealth of 
data at energies around the knee has been accu- 
mulated by a large number of experiments op- 
erated over many years. It is clearly notice- 
able that the data fall into two groups differ- 
ing by their fluxes mostly: CASA-MIA, CASA- 
BLANCA, and DICE (all operated at Dugway, 
Utah) show distinctly lower fluxes than Tibet, 
HEGRA, EAS-TOP, and KASCADE and Tunka 
pS] (not shown in figure 1). This problem has al- 
ready been addressed in [27] but is still not fully 
understood. It may be related to different ob- 
servation techniques (charged particles combined 
either with muons or with Cherenkov light), dif- 
ferences in the details of EAS simulations, or to 
other reasons. On the other hand, it should be 
pointed out that the differences almost vanish, 
if one of the groups is shifted by about 15 % in 
their absolute energy scale, i.e. by an amount 
well within the systematic uncertainties of the 
experiments. The knee energy is found in all ex- 
periments at approximately 3 PeV with the in- 
dex changing from 71 = 2.7 to 72 = 3.1. Only 
Akeno data are different showing different spec- 
tral shapes and a very sharp knee at '--^ 5 PeV. 

It has been realized that the all-particle spec- 
trum is not very discriminative against astrophys- 
ical models of the knee and that a deconvolution 
into different primary particles is required. How- 
ever, this is probably the most difficult task in 
EAS physics, both because of the level of depen- 
dence on hadronic interaction models used in EAS 
simulations and because of the significant (mass 
dependent) fluctuations of EAS observables. A 
large variety of methods is used to infer the pri- 
mary energy and mass [55], most notably the ra- 
tio of electron to muon numbers. At energies 
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Figure 3. CR energy spectrum of light (p-l-He) 
and heavy (all the rest) primaries from EAS-TOP 
and MACRO using TeV muons [.2.9 ). 



higher than approx. 10^^ eV, also direct measure- 
ments of the shower maximum in the atmosphere 
become available by observations of fluorescence 
light with imaging telescopes, such as operated 
by HiRes and the Pierre Auger Collaboration (see 
proceedings to this conference). 

Extensive analyses of both the energy spectrum 
and composition have been performed by EAS- 
TOP and KASCADE. EAS-TOP has analysed its 
data through simultaneous measurements of the 
electromagnetic and muonic shower components. 
These are obtained from the EAS array operated 
at Campo Imperatore on the mountain top 2005 
m.a.s.l. (820 g/cir?) above the underground Gran 
Sasso Laboratories in which the MACRO detec- 
tor has been located under an average depth of 
1200 m rock [29]. The coincident observation of 
the soft charged particles in the surface array and 
the high energy EAS muons {E^ > 1.3 TeV) in 
the underground detector permits - despite large 
fluctuations of the muon number - a reconstruc- 
tion of the CR energy spectrum for "light" and 
"heavy" primaries. The result, depicted in figure 
121 shows that the energy spectrum of the light pri- 
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maries is beginning to diminish at about 5 PeV, 
whilst the heavy component may be signaUng its 
change in the spectrum at least a decade higher 
in energy. 

The results corroborate those of KASCADE 
shown in figure ID KASCADE is located at 
sea-level (ffO m.a.s.l.) in Karlsruhe, Germany, 
and measures the electromagnetic, muonic, and 
hadronic EAS components using a very dense de- 
tector array and a hadronic calorimeter [30j . The 
analysis of the data takes advantage of the ef- 
fect that for given energy, primary Fe-nuclei re- 
sult in more muons and fewer electrons at ground 
as compared to proton primaries. Specifically, in 
the energy range and at the atmospheric depth 
of KASCADE, a Fe-primary yields about 30% 
more muons and almost a factor of two fewer elec- 
trons as compared to a proton primary. The basic 
quantitative procedure of KASCADE for obtain- 
ing the energy and mass of the CRs is a technique 
of unfolding the observed two-dimensional elec- 
tron vs truncated muon number spectrum into 
the energy spectra of primary mass groups |20| . 
The problem can be considered a system of cou- 
pled Fredholm integral equations of the form 
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is another integral equation with the kernel func- 
tion kA = TA • • SA factorizing into three parts. 
Here, ta describes the shower fluctuations, i.e. 
the 2-dim distribution of electron and truncated 
muon number for fixed primary energy and mass, 
describes the trigger efficiency of the exper- 
iment, and SA the reconstruction probabilities, 
i.e. the distribution of and A'^*'' that is recon- 
structed for given true numbers iV* , A^*"''* of elec- 
tron and truncated muon numbers. The proba- 
bilities PA are obtained from CORSIKA simula- 
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Figure 4. Unfolded CR energy spectrum of p, He, 
and C mass-groups from KASCADE. The spectra 
are obtained by using QGSJET and SYBILL for 
the generation of the EAS response matrix pA 

m. 



tions using QGSJET-01 ^ and Sibyll 2.1 [21] 
as high-energy and GHEISHA [31] as low-energy 
hadronic interaction models and a moderate thin- 
ning procedure. Smaller samples of fully sim- 
ulated showers were generated for comparison. 
The simulated data are then fed into the detec- 
tor Monte Carlo programme and the response is 
parameterized as a function of energy and mass. 
Because of the large shower fluctuations, unfold- 
ing of all 26 energy spectra ranging from protons 
to Fe-nuclei is clearly impossible. Therefore, 5 
elements (p. He, C, Si, Fe) were chosen as repre- 
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sentatives for the entire distribution. More mass 
groups do not improve the ^^-uncertainties of the 
unfolding but may result in mutual systematic bi- 
ases of the reconstructed spectra. 

The results of such an unfolding are presented 
in figure [H Shown are the spectra of the p, 
He, and C mass-groups based on the response 
matrices pA obtained from the two interaction 
models. Clearly, there are common features but 
also differences in the energy distributions. In 
each of the distributions a distinct break in the 
spectrum is observed which is increasing towards 
higher energy with increasing primary mass. In 
both cases the He flux is higher than the pro- 
ton flux. This finding may be surprising at first 
sight, but it is already suggested by extrapolat- 
ing the He and proton spectra with their different 
slopes from lower energies towards the knee (see 
figure [2]) • The spectrum of the Si group (see fig- 
ure 14 and 15 in Ref. [5D]) indicates a knee at 
even higher energies. The Fe spectrum (figure [2]) 
shows large differences when performing the un- 
folding either with the QGSJET or SibyU model 
demonstrating that such analyses are limited at 
present mostly by uncertainties of the hadronic 
interaction models. Despite these differences in 
the individual spectra, the all-particle spectra of 
KASCADE (see figure [T]), obtained by summing 
up the energy spectra of all mass groups (p - Fe) 
coincide very nicely for the two interaction mod- 
els. Thus, it can firmly be stated that the knee in 
the all particle spectrum is caused by light (p and 
He) primaries. Obviously, also the mean mass 
composition (e.g. expressed in terms on the mean 
logarithmic mass |32j ) increases above the knee. 

A more detailed investigation [5D] shows that 
the QGSJET model performs reasonably well at 
high energies but exhibits some problems at PeV 
energies. Sibyll, on the other hand, describes the 
data rather well in the knee region but suffers 
from a muon deficit at higher energies. There- 
fore, it suggests a more prominent contribution 
of heavy primaries at high energies. It should 
be emphasized, that this muon deficit of Sibyll 
applies to 0{lGeV) muons only. Muons at en- 
ergies of several 100 GeV, such as observed by 
underground experiments like AMANDA and Ice- 
Cube, seem to be described rather well by Sibyll 
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Figure 5. Interpretation of the CR spectrum in 
terms of different sources [11] . Shown are the in- 
dividual galactic sources (component A and B) 
and the flux expected from extragalactic sources. 
The galactic components are guided by the KAS- 
CADE knee shape as far as the point marked x. 
The dashed line Q is the total if the extended tail 
B of the galactic flux is omitted. 



[33] . Very recently, a new interaction model, 
called EPOS has been released [31]. Most impor- 
tantly, it provides a better description of baryon- 
antibaryon production at high energies. A pre- 
liminary analysis shows that the muon num- 
ber increases more rapidly with energy than in 
QGSJET or Sibyll with the muon density be- 
ing about 40 % higher at 10^^ eV compared to 
QGSJET-01 calculations. It will be interesting 
to repeat the unfolding of the CR energy spec- 
tra employing this model to verify whether the 
present deficiencies of the interaction models will 
be resolved. 

The unfolded KASCADE energy spectra can 
directly be compared to phenomenological calcu- 
lations of astro- and particle physics related mod- 
els or can be used to infer information about the 
CR sources. An example is shown in figure (5) 
taken from Ref. [11]. It is concluded that the 
data provide support for the supernova picture 
of CR origin, i.e. the distinct knee near 3 PeV 
would be related to emission by the free expan- 
sion phase of SNRs. However, a question arises 
about how to fill the gap from the iron knee at 
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about 10^^ eV to the ankle at ~ 5 ■ 10^^ eV. 
These CRs may originate from SN type II ex- 
plosions into dense stellar winds where the inter- 
action generates much stronger magnetic fields. 
This may result in rigidities up to at least 10^^ V 
(component 'B' in figure[5]), especially from a few 
abnormally high speed/low mass ejections [11] . 

A very important question is whether the 
present data allow to distinguish a knee of con- 
stant rigidity {E/Z) from that of constant energy 
per nucleon (E/A), such as is predicted by parti- 
cle physics interpretations of the knee or by the 
cannonball model. Unfortunately, Z/A changes 
only from 0.5 in case of He to 0.46 for Fe nu- 
clei. Hence, the question about the rigidity de- 
pendence needs to be answered basically by com- 
paring the energy spectra of p and He primaries. 
Ironically, these are the two primaries which are 
most strongly affected by EAS fluctuations, so 
that their energy resolutions are deteriorated sub- 
stantially. In fact, overlaying the p and He spec- 
tra of figure [4] using E/Z and E/A abscissas does 
not give a clear answer; Sibyll exhibits a slight 
preference for charge scaling and QGSJET for 
mass scaling. It is hoped, that the situation will 
improve somewhat with better models becoming 
available. Improving on the data side seems more 
difficult because of two reasons: statistical errors 
are already much smaller than systematical ones 
and (presently not yet included) data from larger 
zenith angles are subject to even stronger EAS 
fluctuations. 

3. THE SECOND KNEE AND ANKLE: 
TRANSITION FROM GALACTIC TO 
EXTRAGALACTIC COSMIC RAYS 

Besides the prominent knee in the all-particle 
spectrum, additional structures are observed at 
about 10^^ eV and ^ 3 • 10^^ eV, known as the 
second knee and the ankle, respectively (see Figs. 
[Hand [5]). The ankle has been reported convinc- 
ingly by a number of experiments, but there is 
still no consensus about the existence of a sec- 
ond knee. This is because of both the weakness 
of the structure making it difficult to detect and 
because of only few experimental data, most of 
which are either at their upper or lower limit of 
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Figure 6. All-particle CR energy spectra from 
Yakutsk [35], Haverah Park [37], Fly's Eye [33], 
HiRes II (mono) [33], and Akeno [35] , 



detectable energies. A blow-up of the data be- 
tween 10^^ eV and 10^^ eV is shown in figure 
[Hi It includes measurements by Akeno |3S] , Fly's 
Eye (stereo) [36], Haverah Park [37], Yakutsk 
and HiRes II (mono) [39] . Akeno has provided the 
first hint of a change in the index of the power-law 
energy spectrum around 6 • 10^^ eV. The steep- 
ening of the spectrum was confirmed by Haverah 
Park and is indicated also in the Fly's Eye and 
more recent HiRes data. A recent re-analysis of 
the Yakutsk 1974-2004 data agrees well with the 
Akeno data providing additional support for the 
existence of a second knee at about (6 ± 2) • 10^^ 
eV. The ankle at ~ 3 • 10^^ eV was first ob- 
served by Haverah Park, Akeno, and Yakutsk and 
is traditionally explained in terms of the transi- 
tion from galactic to EGCRs. The key point here 
is that one expects the galactic magnetic field to 
lose its efficiency at about this energy as the gyro- 
radius of a particle at charge Z ma. /^-Gauss field, 
~ 1 kpcZ'~^_B~Q, becomes comparable to the 
thickness of the galactic disk. It then becomes 
natural to think of hard EGCRs starting to pen- 
etrate into the galaxy and dominating the fiux at 
higher energies (see figure [5] for illustration). 
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Figure ini also provides an intuitive explanation 
for the second knee: it would primarily be caused 
by the break of the galactic Fe component. As 
Hillas pointed out [11] . an extra component 'B' 
would be needed in order to make up the well- 
measured total CR flux at several 10^^ eV for 
which he considered SNae Type II explosions into 
dense stellar winds (see chapter 2.2). Naively, the 
second knee in this picture is expected at Epc — 
26 X 3 • 10^^ ~ 8 • 10^^ eV or even lower if the knee 
is composed of p and He primaries as suggested 
by figure ID This is almost a factor of 10 lower 
than reported by Akeno and others. A scaling of 
the knee position with E/A would bring the Fe- 
knee up to approx. 2 • 10^^ eV, but still too low 
to fit the classical picture. 

Ignoring this puzzle for a moment, also char- 
acteristic changes of the CR composition are ex- 
pected in this traditional picture of the knees and 
ankle. Up to the knee, the composition would fol- 
low the standard source composition dominated 
by p and He primaries. Between the first and 
second knee the composition would change to 
become iron dominated, and above the ankle it 
would be dominated by extragalactic protons. 

However, the 'folklore' about the second knee 
and ankle and its related transition from galac- 
tic to EGCRs is not free of dispute and has 
received much attention recently. Back in the 
80s, Berezinsky and collaborators have pointed 
out an inevitable feature of the 10^^-10^^ eV 
EGCR spectrum: if EGCRs consist of protons 
mostly, they would suffer - besides the GZK ef- 
fect - from energy losses associated with the pro- 
duction of e+e" pairs in the CMB photon field 
[40] . This would result in a modulation of the 
all-particle energy spectrum to what is called a 
"pair-production dip" between 1 ■ 10^* - 4 • IQ-'^^ 
eV. In such a way, the turn over from the left- 
to the right hand side of the 'dip' would mimic 
the ankle. Moreover, since the Bethe-Heitler pair 
production works effectively only for protons [41j , 
the ankle can then be interpreted as a signature of 
a pure proton EGCR component and the galactic- 
extragalactic transition must occur at much lower 
energies than in the traditional picture, possibly 
around the second knee. 

How can the two models be discriminated? The 
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Figure 7. Comparison of the mean depth of 
shower maximum, Xmax, predicted by the dip- 
and ankle-model of Ref. and the Cannonball 
of Ref. [T^ with data from HiRes 03], [Mj, Fly's 
Eye [36] , and Yakutsk [45] (see text for details) . 



most critical observation is provided by a mea- 
surement of the chemical composition in the en- 
ergy range around 10^® eV: In the dip model a 
strong dominance of protons, and in the ankle 
model a strong dominance of iron nuclei is ex- 
pected. A recent confrontation of the two models 
to existing data has been performed by AUard 
et al. [42]. The authors conclude that the all- 
particle energy spectrum is reproduced equally 
well by the two models. However, based on a com- 
parison of the mean mass composition, analyzed 
in terms of the mean depth of the shower max- 
imum, ATmax, they favour the traditional model. 
Figure [7| compares the ATmax data of various ex- 
periments with the dip- and ankle- model of [42j . 
Here, only CORSIKA / QGSJET-01 simulations 
are shown, because QGSJET-01 is the interac- 
tion model providing the most consistent descrip- 
tion of experimental data in this energy region. 
Clearly, this direct comparison with QGSJET-01 
does not seem to give preference to any of the two 
models. Also shown are predictions of the can- 
nonball model for two choices of penetrability of 
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EGCRs into the Galaxy [TU]. It should be noted, 
that in the latter case, X^i,^{E) is constructed by 
a simplified model described in Ref. HHj instead 
of using full EAS simulations. Evidently, better 
data are required before definite conclusions can 
be drawn about the transition from Galactic to 
extra-galactic CRs. 

4. ANISOTROPIES 

Another key observation in cosmic ray astro- 
physics is the directional distribution of the parti- 
cles. That distribution will depend on any galac- 
tic magnetic fields and hence will be energy (rigid- 
ity) dependent. However, with very limited ex- 
ceptions, which are not individually statistically 
significant, there is no observed deviation from 
isotropy above the knee of the energy spectrum, 
and any anisotropics at lower energies are them- 
selves very small 47 48 . Probably, the most com- 
prehensive data at energies from a few to several 
hundred TeV have been obtained by the Tibet 
AS7 experiment. Besides revealing fine details of 
known anisotropics, the data support the picture 
of corotation of low energy CRs with the local 
Galactic magnetic environment and they may in- 
dicate an anisotropy around the Cygnus region 
|49j . However, a contamination of TeV 7's in the 
data sample cannot be excluded at present. 

A non-uniform distribution in of the arrival di- 
rections, suggestive of a source direction, in the 
energy range 10^^ ° - 10^^ '* eV has been reported 
by the AGASA [50] and similarly by the SUGAR 
collaboration f5T\. However, neither of those ob- 
servations on their own is clearly statistically sig- 
nificant. Moreover, the Pierre Auger Collabora- 
tion has also started to analyse the galactic centre 
region. These results, obtained with much larger 
exposure than of AGASA and SUGAR, do not 
support that finding and instead provide an up- 
per bound on a point-like flux of CRs from the 
Galactic Centre. Even in absence of CR point 
sources, such data may be regarded as the pos- 
sible beginning of a new era in cosmic ray astro- 
physics in which we can begin directional cosmic 
ray astronomy. The possibility of having a source 
to observe may indeed open up new frontiers for 
the Pierre Auger Observatory |53I54| . 
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Figure 8. Rayleigh amplitudes A vs primary en- 
ergy from different experiments. The data of 
KASCADE gH] (bold line) represent upper lim- 
its (95 %). The thin lines show expectations from 
the galactic CR diffusion model of Ref. |55j . 



As already pointed out, the low level of CR 
anisotropy even at energies above the knee is con- 
sidered the most serious challenge to the stan- 
dard model of the origin of galactic CRs from 
diffuse shock acceleration [11] . Figure [8] shows 
a collection of data expressed in terms of the 
Rayleigh amplitudes A. The thin lines represent 
a CR diffusion model [55] predicting anisotropics 
on a scale of 10"'' to 10~^ depending on particle 
energy and strength and structure of the galac- 
tic magnetic field. However, the model fails to 
describe the all-particle spectrum considerably. 
Assuming a simple rigidity model of TcsciE) oc 
^-0.6^ Hillas estimates anisotropies at a level of 
5%, 16%, and 180% at 1.5 • lO" eV, lO^^ eV, 
and 1.5 • 10^^ eV, respectively [11]. In case of a 
E^^/"^ scaling, the values would go down to 0.6 %, 
1.1 %, and 3.7% which is still in contradiction to 
the experimental data of figure HI 

As already mentioned, the CB-model [TD] pre- 
dicts much lower levels of anisotropies than mod- 
els in which CRs diffuse away from the central 
realms of the Galaxy, where most SN explosions 
take place. A CB, on the contrary, is consid- 
ered a continuous source of CRs along its trajec- 
tory from the galactic disk into the galactic halo. 
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Along its trajectory, the source intensity depends 
on the local and previously traversed ISM den- 
sity. Thus, the source of CRs is very diffuse and 
the directional anisotropy of CRs at the Earths 
location is expected to be very small and to vary 
little with energy. 

5. SUMMARY AND OUTLOOK 

Diffusive shock acceleration in SNRs is consid- 
ered a viable mechanism for accelerating cosmic 
rays and it naturally leads to a power-law spec- 
trum in rigidity. However, many fundamental 
questions related to the assumption of SNRs be- 
ing the sources of galactic cosmic rays arc still 
open. These questions include, amongst others, 
the absence of TeV 7-radiation from a large frac- 
tion of SNRs, the origin of the knee in the cosmic 
ray spectrum, the low level of global anisotropics 
in their arrival direction, the transition of galac- 
tic to extragalactic CRs and its related question 
about the existence of a second knee and about 
the origin of the ankle. 

The deconvolution of the all-particle CR spec- 
trum into energy spectra of individual mass- 
groups by current KASCADE data 20 has ad- 
vanced the field quite a lot. Such kind of data con- 
tain much more information than the all-particle 
spectrum and the mean mass of CRs (expressed 
mostly by (ATmax) and (In A)) alone. However, 
there remain large uncertainties, which still allow 
alternative interpretations. Most prominently, a 
definite answer about an EjZ (rigidity) or EjA 
scaling of the knee position cannot be given at 
present. However, there is still some room for 
improving the data quality and, despite enor- 
mous progress already made, there are also bet- 
ter hadronic interaction models being developed 
which arc hoped to eliminate the still existing 
shortcomings of the present models such as Sibyll 
2.1 or QGSJETOl [S]. 

At energies above 10^^ eV data become very 
sparse and we are far from understanding the 
transition from galactic to extragalactic CRs. Al- 
though the ankle in the CR spectrum at about 
5 • 10^* eV is often interpreted as the signature 
of the transition from a steeply falling galactic 
CR-spectrum to a slightly harder extragalactic 



spectrum, alternative explanations are possible. 
Sometimes the second knee at about 10^^'^ cV is 
considered as indication for the transition to ex- 
tragalactic CRs, but this explanation would re- 
quire fine-tuning of the injection spectra of the 
different galactic and extragalactic sources. Two 
particular models were discussed in detail, the 
dip- 41j and the ankle-model |Tl]. Current data 
on (Amax) do not allow to exclude any of the two 
models. The transition from galactic to extra- 
galactic CRs occurs in the energy region of the 
second knee and is distinctly seen only if iron and 
proton spectra are measured separately. 

In conclusion, the fundamental question about 
the origin of high energy CRs below the GZK en- 
ergy remains far from being answered. As a con- 
sequence, the interest in studying CRs from about 
10^^ to 10^^ eV with high quality state of the art 
EAS detectors has grown worldwide and several 
new experiments are being prepared or planned 
for. These include KASCADE-Grande (already 
in operation) [56] as well as low-energy extensions 
of Auger by High Elevation Auger Telescopes 
(HEAT) and an infill array with extra muon de- 
tectors, as well as the Telescope Array (TA) and 
its low-energy extension TALE 57J. These detec- 
tors can reliably solve the problem of measuring 
the energy spectrum and mass composition in the 
transition region and complement the measure- 
ments performed at the highest energies by the 
Pierre Auger Observatory. 
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